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ELECTRIC WINPTNG DTS PTACEMENT PETECTTQN MP/rrrnn Aism 
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FIELD OF THE INVENTION 



5 



The present invention relates generally to conductive windings used in electrical 



components. More particularly, the present invention relates to conductive windings found in 
transformers and electrical motors and generators. 



one potential level to another. 

Transformers are one of the most essential elements of an electric power system. They 

are widely used in electric networks because the generation, transmission, and distribution of 

power require different voltage levels. For example, in a large power utility, the number of 
15 transformers can exceed 1000. They link electric loads to power supplies through an 

interconnected power network and satisfy the requirem^ts of both parts. The use of 

transformers helps to reduce losses in an AC power system. 

Transformers are one of the most expensive pieces of equipment in a power grid. A 

large power transformer in a 500kV system may cost in the order of a million dollars. In 
20 1996, the total sale of power transformers in the United States was over US$ 1 1 million. In 

temis of numbers, more than 8000 transformers, of at least 2MVA and up to more than 100 

MVA, were sold in the same year all over the world. 

Since transformers play a vital role in the operation of AC or HVDC electric networks, 

it is important to ensure that they are operating efficiently and reliably. Furthermore, the 
25 failure of a transformer can be potentially devastating to personnel safety and to the 

environment. Transformers may fail explosively causing great personal injuries to people 

around them and damage to the surrounding equipment. The environment can be adversely 



affected by the leakage of oil. A failure also has a large economic impact due to its high cost 
of replacement and repair and lost revenue while it is out of service. During its outage, the 
30 customers of a utility can be greatly inconvenienced, which could result in a loss of goodwill 
towards the utility. 
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In an AC power system, power transformers are used to convert electric power from 
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Transformers are well known in the art» and can include dry type transformers and 
fluid-filled transformers, with different core configurations. Figure 1 illustrates an iron-core 
transformer which is typically encased within a tank. Generally, the iron-core transformer 
consists of one or more sets of windings that are coupled through a common magnetic core. 
5 The main components of transformer 10 shown in Figure 1 include a laminated iron core 12, 
and a set of windings 14a and 14b, where winding 14a can be a high voltage winding while 
winding 14b can be a low voltage winding. Each set of windings normally consist of two or 
three windings: a primary winding, a secondary winding, and sometimes a tertiary winding, 
mounted concentrically on top of each other. As shown in Figure 1, windings 14a and 14b 

10 are wrapped around the legs of the rectangular window iron core 12. The windings 14 are 
usually made of copper or aluminum, and can be shaped as either wires or sheets with an 
insulating layer 16 separating successive layers of the windings. The windings 14 are 
clamped to the tank structure to maintain their physical positioning. Core clamps 18 are used 
to fix the core to the tank structure. Ceramic bushings are used to isolate the windings from 

15 grounded structures of the transformer such as the oil tank. Mineral oil is typically used as 
insulation medium, and for cooling the transformer. In operation, a time-varying flux created 
by one winding induces a voltage in the other winding. Figure 1 represents one known 
configuration of a core-type transformer, and those of skill in the art will understand that there 
are a number of laiown configurations. 

20 Inevitably in service transformer failures occur due to either internal failures and/or 

external failures. As the context of the present invention is related to internal failures, 
external failures will not be discussed. 

Internal failures are faults that occur inside the tank, such as short circuits between 
windings, short circuits between windings and iron, short circuits between turns, insulation 

25 deterioration, loss of winding clamping, partial discharges and winding resonance. One of the 
primary causes of internal failures is winding movement, or displacement, which can lead to 
insulation deterioration, and winding collapse. When the insulation deteriorates, windings 
can electrically contact each other, resulting in high current flow, partial discharges, and very 
severe and expensive faults. 

30 An external short circuit fault, for example, created by lightning strike or ground £ault, 

is the most likely factor to cause winding movement. Assuming a copper winding is wound 
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on a ferromagnetic core and the transformer is in service, the currents carried by the windings 
will produce the predominant flux in the axial direction. The interaction of the current in the 
coils in the circumferential direction and the axial flux field will therefore produce a radial 
outward electro-magnetic force. Figure 2 illustrates possible movement of a winding 20 
5 relative to a core 22 of a typical iron-core transformer. A current I flowing through the 

winding 20 in the direction as indicated by the arrows^ can suffer from movement in the axial 
direction as shown by arrows 24, and outward winding movement in the radial direction as 
shown by arrows 26. There is, however, a leakage flux around individual turns and near the 
ends of the windings which has a component in the radial direction. This radial flux induces 

10 an axial inward electro-magnetic force. Under normal operation of power transformer, the 

windings are designed to withstand the mechanical pressure described above. However, when 
short circuits happen, the electro-magnetic forces induced in the windings are increased 
dramatically and threaten the insulation layers severely. For example, if a transformer's 
leakage impedance is 10%, its short circuit current will be 10 times the rated current and the 

15 mechanical stress will be roughly 100 times of the normal stress under the rated load current. 

In the event of a short circuit situation, the winding will stretch out in the radial 
direction and compress in the axial direction. The radial forces in a lOMVA transformer can 
exceed 100,0001bs. Such huge electro-magnetic forces will inevitably loosen the winding 
clamp structure, which is the main mechanical support of the winding and cause the distortion 

20 or movement of the winding. 

hirush currents, due to the powering up of the transformer circuit, and vibration forces 
increase the electro-magnetic forces in the same way as short circuits. The frequent 
fluctuation in generation or load will also put burdens to the electrical and mechanical 
strength of the windings, accelerate the loosening of the winding clamps, and eventually cause 

25 the winding to move or even break down. Furthermore, winding bulges or sharp edges of the 
coil can cut through insulation and cause short circuits between turns, which can be 
exacerbated as the winding moves. 

Since a large number of failures are due to the windings, techniques for diagnosing the 
health of a transformer have been proposed, since it is useful to be able to assess the health of 

30 a winding to enable prediction of remaining life, capacity limits, and preventive maintenance. 
A short summary of some of these methods follows. 
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There are several known methods to assess the health of a wmding, including 
detection of insulation degradation. By way of example, insulation degradation in 
transformers can be monitored by methods that detect partial discharges. Partial discharge 
detection methods include Dissolved Gas Analysis and Tan-Delta techniques. It is notable 
5 that Dissolved Cras Analysis methods detect the chemical reaction products of insulation 
degradation, and therefore partial discharge must have occurred before these techniques are 
useful. The Tan-Delta similarly detects damage that has already occurred. 

The winding ratio test measures the numbers of turns of both primary and secondary 
windings and calculates the ratio between them. By comparing the measured winding ratio 

10 with the ratio of rated primary and secondary voltages as shown on the nameplate of the 

transformer, shorted turns or open winding faults may be detected. However, an outage and 
isolation of the transformer is required for the pmpose of measurements. 

The winding resistance test, which is similar to the winding ratio test, except that it 
measures the winding resistance rather than the number of winding turns. Additionally, a very 

15 precise ohmmeter is needed, wliich will assure the accuracy of a fraction of an ohm. The 
measured resistance will be compared with the previous measurement referred to the same 
temperature. Measurements are conducted for different phases and different tap-changer 
positions. This method detects the condition of the winding conductor directly. However, it 
requires a transformer outage and is usually performed in the factory or a laboratory. 

20 The 60 Hz transformer impedance test measures the input voltage, input current and 

input power while shorting the low voltage winding of the transformer. The before and after 
results are then compared with the before and after short circuit. This test is insensitive to 
small winding movements. 

The leakage reactance measiurement (LRM), which can be achieved with the same test 

25 set-up as the short circuit impedance measurement, is based on the increased leakage 

reactance resulting firom the radial outward force on the outer winding and the radial inward 
force on the imier winding induced by short circuits. 

Another way to assess the healdi of a winding is to detect winding displacement. 
Detection of winding displacement is advantageous because such methods permit detection of 

30 deterioration prior to actual winding damage. Known winding displacement detection 

methods include detection of increases of the audible noise, visual inspection; short-circuit 
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impedance measurement, vibration analysis, low voltage impulse, and swept Frequency 
Response Analysis ('FRA') techniques. 

The FRA technique compares die input admittance function Y(co) of the displaced 
winding with the Y(a)) of an equivalent healthy winding. This technique is well known in the 
5 art, and looks at the transformer as a lumped impedance and measures the admittance function 
Y(w) as the ratio of input current to applied voltage for a range of high frequencies (e.g., 1 to 
10 MHz). The idea of the FRA test is to compare the "signature" of the transformer, as 
defined by the shape of the Y(ca) function, as it changes with aging. 

Figure 3 shows measurement plots made using the FRA technique of a healthy 

10 winding 20 versus that of a displaced winding 22. The admittance Y(co) in Siemens is plotted 
on the vertical axis while frequency in MHz is plotted on the horizontal axis. Two main 
problems exist with this technique. First, the frequency response differences between the two 
curves only begin to appear at frequencies beyond 2-4 MHz. Unfortunately, online 
measurements are corrupted due to electrical noise interference, which becomes problematic 

15 at these frequencies. Second, the observed response distortion is difficult to directly correlate 
any specific problem with the displaced winding, and in particular, to the amount of winding 
displacement. 

The main problem with the FRA technique is that the Y(co) function is highly 
oscillatory because it includes all the winding electrical parameters: resistance, inductance, 

20 and capacitance and, in addition, is very sensitive to electrical noise and to numerical noise in 
the processing of the signals. As a result, FRA tests are normally performed off-line with the 
equipment removed from service and in a very controlled test environment. Taking a power 
transformer out of service is very expensive and sometimes not even feasible because of 
service continuity constraints. 

25 It is, therefore, desirable to provide a method and system for assessing the health of 

electrical windings in transformers without having to remove the transformer unit from 
service, while providing effective determination of winding displacement and/or faults in the 
winding. 
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SUMMARY OF THE INVENTION 

It is an object of the present invention to obviate or mitigate at least one disadvantage 
of prievious electrical winding health assessment techniques. Jn particular, it is an object of the 
present invention to provide a method and apparatus for effectively assessing Ihe health of an 
5 electrical winding with high sensitivity and without having to disconnect or remove the unit 
containing the winding from its environment. 

In a first aspect, the present invention provides a diagnostic device for generating a 
characteristic impedance for an electrical winding having an input terminal and an output 
terminal. The diagnostic device including a signal generator, sensing means, processing 

10 means and storage means. The signal generator applies an electrical signal having a frequency 
component to the input terminal of the electrical winding. The sensing means detects a 
magnitude and a phase of an output electrical signal at the output terminal of the electrical 
winding. The sensing means then converts the magnitude and the phase of ttie output 
electrical signal into digital signals. The processing means sets parameters of the electrical 

15 signal and receives the digital signals. The processing means then calculates the characteristic 
impedance with the digital signals and the parameters of the electrical signals corresponding 
to one frequency, based on a transmission line model of the electrical winding. The storage 
means stores the digital signals and the parameters of the electrical signal. 

In embodiments of the present aspect, the signal generator can include one of a 

20 function generator and a network analyzer for generating the electrical signal at the one 
frequency, the one frequency can be a minimum frequency selected corresponding to a 
maximum length of the electrical winding, or the one frequency is at least about 500kHz. The 
sensing means can include a high speed digital data recorder, and the storage means can 
include a memory device. 

25 According to another embodiment of the present aspect, the processing means 

includes a sequence generator and a calculation engine. The sequence generator sets the 
parameters of the electrical signal for generation by the signal generator, where the 
parameters can including input voltage, input current and the one frequency. The calculation 
engine receives the input voltage, the input current, the one frequency, and the digital signals, 

30 and is programmed with the transmission line model of the electrical winding for calculating 
the characteristic impedance. The processing means can further include a controller core, a 
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fuzzy logic engine and a graphic engine. The controller core can receive Hie characteristic 
impedance from the calculation engine, for comparing the characteristic impedance with a 
base characteristic impedance to provide a corresponding difference value. The fuzzy logic 
engine can receive the corresponding difference value and apply pre-programmed fiizzy logic 
5 to provide a corresponding text message. The graphic engine can receive the characteristic 
impedance and the base characteristic impedance for generating graphical plot information. 

In further aspects of the present embodiment, the controller calculates an approximate 
displacement value of the electrical winding corresponding to the difference value, and the 
user interface displays the corresponding text message and displays the graphical plot 
10 information as a graph of impedance versus frequency. 

In yet another embodiment of the present aspect, the transmission line model of the 
electrical winding is expressed by 

15 

for representing the circuit model of the electrical winding shown in Figure 6. The 
transmission line model of the electrical winding can be further expressed by 

(V,+ZMe-^^ =V^+ZJf; 

where each expression represents a loop in the circuit model of the electrical winding shown 
in Figure 8. 

In further embodiments of the present aspect, the signal generator can include a 
25 recurrent surge generator for generating the electrical signal, where the electrical signal 

includes a train of pulses, such as a train of square pulses. The processing means can include a 
calculation engine for executing a Fourier transform algorithm to decompose the electrical 
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signal and the digital signals into frequency components. The electrical winding can include a 
transfonner winding housed in a tank, and the transformer can be on-line. 

In a second aspect, the present invention provides a method for determining a 
characteristic impedance of an electrical winding. The method includes the steps of applying 
6 an input signal having a frequency component to a first terminal of the electrical winding; 
measuring an output signal at a second terminal of the electrical winding; storing the input 
signal data and the output signal data; and calculating the characteristic impedance based on a 
transmission line model of the electrical winding with the input signal data and the output 
signal data corresponding to one frequency. 

10 According to embodiments of the present aspect, the step of applying can include 

generating an analog signal having predetermined voltage and current values at the one 
frequency, the one frequency can be a minimum frequency selected corresponding to a 
maximum length of the electrical winding, or the one frequency can be at least about SOOkHz. 
In another embodiment of the present aspect, the step of applying includes generating 

15 a pulse train having predetermined voltage and current values, and the step of calculating 
includes executing a Fourier transform algorithm for decomposing the input signal and the 
output signal into frequency components, where the frequency components include the one 
frequency. The one frequency can be at least about SOOkHz. 

In yet another embodiment of the present aspect, the steps of applying, measuring and 

20 storing are repeated for a plurality of distinct frequencies before the step of calculating. The 
step of calculating can include calculating the characteristic impedance of the electrical 
winding at each of the plurality of distinct frequencies. The electrical winding can include a 
transformer winding housed in a tank, and the transformer can be on-line. 

In a third aspect, the present invention provides a method for assessing a condition of 

25 an electrical winding. The method includes a) obtaining a base characteristic impedance of the 
electrical winding at a first time; b) storing the base characteristic impedance corresponding to 
the first time; c) obtaining a current characteristic impedance of the electrical winding at a 
second time after the first time; d) storing the current characteristic impedance corresponding 
to the second time; e) calculating a difference value between the current characteristic 

30 impedance and the base characteristic impedance; and £) calculating an approximate winding 
displacement from the difference value. 
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16 



20 



According to embodiments of the present aspect, the mefliod further includes a step of 
calculating an approximate fault impedance of the electrical winding at the second time, and 
the step of obtaining the base characteristic impedance includes calculating the base 
characteristic impedance from the expression 



Zc is the base characteristic impedance, 

ju is a magnetic permeability of the electrical winding material 

€ is a dielectric constant of insulation material 

N is a total number of the turns of the electrical winding 

h is an axial length of the electrical winding 

h is an insulating distance between the electrical winding and a core of a transformer 

d is an insulating distance between the electrical winding and a tank of the transformer 
Vois the velocity of light in a vacuum, 300m/us. 

In yet another embodiment of the present aspect, the step of obtaining the base 
characteristic impedance includes applying an input signal having a frequency component to a 
first terminal of the electrical winding; measuring an output signal at a second terminal of the 
electrical winding; storing the input signal data and the output signal data; and calculating the 
base characteristic impedance from a transmission line model of the electrical winding with 
the input signal data and the output signal data corresponding to one frequency. The steps of 
applying, measuring, storing and calculating can be repeated for a plurality of frequencies. 

hi a further aspect of the present embodiment, the step of obtaining the current 
characteristic impedance can include the steps of applying the input signal having the 
frequency component to the first terminal of the electrical winding; measuring a second 
output signal at the second terminal of the electrical winding; storing the input signal data and 
the second output signal data; and calculating the current characteristic impedance from the 
transmission line model of the electrical winding with the input signal data and the second 
output signal data corresponding to the one frequency. The steps of applying, measuring, 
storing and calculating can be repeated for the plurality of frequencies. 




M N b 
€ h 1 + b/d 



where: 



-9- 



wo 20D5/091004 



PCT/CA2005/000438 



Other aspects and features of the present invention will become apparent to those 
ordinarily skilled in the art upon review of the following description of specific embodiments 
of the invention in conjunction with the accompanying figures. 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

Embodiments of the present invention will now be described, by way of example only, 
with reference to the attached Figures, wherein: 

Fig. 1 is an illustration of an iron core transformer; 
Fig. 2 is an illustration of winding movement in an iron core transformer; 
10 Fig. 3 is a plot of measxirements of windings with and without winding 

movement, obtained using the FRA technique; 

Fig. 4 is a circuit diagram of a transmission line; 
Fig. 5 is a circuit diagram of a frequency-dependent line model; 
Fig. 6 is a circuit diagram showing connections of a winding under 
15 investigation; 

Fig. 7 is a circuit diagram of a winding model; 
Fig. 8 is a circuit diagram of a faulted winding model; 
Fig. 9 is Matlab script for solving Zc equations; 
Fig. 10 is a plot of Zc for different winding movements; 
20 Fig. 11 is a plot of Zc for different winding faults; 

Fig. 12 is a block diagram of a winding diagnostic device according to an 
embodiment of the present invention; 

Fig. 13 is a block diagram of the functional blocks of the processing means 
shown in Figure 12; 

25 Fig. 14 is a flow chart showing a method of calculating Zc for a winding; and. 

Fig. 15 is a flow chart showing a method of diagnosing the health of a 
winding. 

DETAILED DESCRIPTION 

30 Generally, the present invention provides a method and system for obtaining the 

characteristic impedance of an electrical winding by measuring the input and output voltages 
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and the currents within a winding across a lower firequency range than the FRA technique, 
and applying transmission line properties to model the winding. The characteristic impedance 
(Zc) is directly proportional to the capacitance of the winding, and is independent of external 
circuits. Thus any changes to Zc will reflect movements of the winding that would affect the 
5 capacitance. Because Zc has a smooth and robust monotonical relationship with frequency of 
the applied signal, the resulting curve is exponential in shape. A relative comparison between 
Zc curves will provide a clear indication of the overall axial or radial winding movements, 
which can be used to assess the ovemll health of the winding. 

The following embodiments of the present invention utilizes the wave propagation 

10 property and frequency dependency of a transmission line model to obtain a characteristic 
impedance Zc of an electrical winding, be it a winding of a transformer, electrical motor, of 
electrical generator. This is due to the fact that the windings of a transformer, electrical motor, 
or electrical generator are sufficiently long to be treated as transmission lines. 

A traditional long-wire device, such as an overhead transmission line, is considered to 

15 be "electrically long" based on its physical lengdi, and frequency of the 

signal being carried by the line. A wire with a physical length / that accommodates several 
wavelengths / of a signal travelling at velocity a can be considered a propagation channel and 
can be modelled using transmission line equations. 

The presently described technique for testing of transformer windings is based on the 

20 realization that a transformer winding is a long-wire device, and at high frequencies, becomes 
a wave propagation channel. As an example, the winding of a typical power transformer can 
be about 2000 meters long. Signals propagate at the speed of light in the medium, which will 
be typically from 200 to 300 km per second. Given a signal of 1 MHz has a period of 1 
microsecond, 2000 meters of winding will have about ten wavelengths of the signal, which 

26 are sufficient to characterize the winding as a propagation channel and allow the use of wave 
propagation functions for its description. 

Following is a discussion of the derivation of the expressions for modeling a 
transformer winding, which is subsequently used for calculating the characteristic impedance 
Zc of a winding. 
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Figure 4 shows a typical transmission line, for the purpose of illustrating the following 
derivation of its frequency-dependent model. The following wave equations for the 
transmission line of Figure 4 are as follows: 

X ^ Jk^ ^ ^ bk^ (1) 

5 

1 1 

where: 

.Vf=ZJf (3) 

V„ = -ZJ, (4) 
10 = V characteristic impedance of the line (5) 

is per imit length propagation constant of the line (6) 

Z' is the per unit length impedance of the line, and 
y is the per unit length admittance of the line. 

15 By combining equations (1) and (2), a forward perturbation function can be found as 

given in equation (7) under the condition + Z^/^ = at x=0. 



20 



30 



where: 



V,-^ZJ^=iV,+ZJ,)e-^ (7) 

Therefore, at node m, equation (7) becomes: 



V„-Z,Il, = (V, + ZJ,)e-r^ (8) 



25 where: 



I is the length of tlie transmission line 



(9) 
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An equivalent circuit as shown in Figure 5 follows from equations (8) and (9). 
Ef^f^ and E^^^ are history functions and are defined as: 
5 = (F, + ZJ,)e-'' (10) 

£,, = (.V„-ZJ„,)e-'' (11) 

From Figure 5 and the above equations, a complicated distributed frequency- 
10 dependent transmission line can be modeled as two separate lumped circuits. Each circuit 

includes the characteristic impedance of the transmission line, the wave propagation constant, 
and a history function. For a given transmission line, its characteristic impedance and 
propagation constant are fixed, since they are determined by the geometry and composition of 
the transmission line. Typical Zc of a transmission line is a smooth almost exponential curve, 
1 5 which decreases with frequency. 

Following is a discussion of the derivation of the expressions for modeling a non- 
faulted transformer winding. The winding of a transformer behaves as transmission line if the 
frequency of the exciting voltage is sufficiently high. An individual winding is modeled by a 
black box which has a specific structure. The previously presented frequency-dependent 
20 model is used for the structure of the black box. The following discussion details the 

derivation of the parameters of the black box given that the input voltage, input current, and 
output voltage and output current are known or measured. The unknown parameters of the 
model are the characteristic impedance Zc and propagation constant of the equivalent line 
model. 

25 Reference is made to Figure 6, which illustrates the model of the winding under 

investigation, as it would be physically connected. The transmission line model for the 
winding as connected in Figure 6 is shown in Figure 7. Applying Kirchoffs voltage law 
around the voltage loop of the right-hand side circuit of Figure 7 gives: 

30 E„,=iZ,+R)I,„ (12) 

where: 
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= I^AO^ is the measured output current 

Substituting the definition of E^^^ from Equation (12) yields: 

V,e'' ^ZJ,e' ^ZJ^^RI^^ (13) 

5 

where: 

- E^Z, 0^ is the measured input voltage 
-If^/LO is the measured input current 

10 Equation (13) contains two complex unknowns, thus another equation must be 

obtained. The second equation comes from applying Kirchoffs voltage law around the 
voltage loop of the left-hand side circuit as given by Equation (14): 



15 



Vk^E^H'^ZJ, (14) 

Using the definition of Ef^^^ from equation (13) and noting the direction of yields: 
V, -ZJ, =RI^e-r -ZJ,„e-' (15) 

Equations (13) and (1 S) can now be solved numerically to obtain the values for the 
20 characteristic impedance Zc and propagation constant / • Note fliat the propagation constant 
r obtained by this method is the actual value not the per-unit length value, i.e. y = y^l . While 

from formula (16), r can be obtained, which is the travelling time of wave propagation 
corresponding to certain co . Formula (17) enables identification of the equivalent length of 
the winding, i.e 7^, which is in between the physical length of the winding and the height of 

25 the winding. Vq is the velocity of wave propagation, which is 3e'*'^m/s in air. 

a> * T = imaginary _ part(j) (16) 
/«=Vo*ir (17) 

30 From the above, it can be seen that by taldng measurements of currents and voltages 

on the two ends of the winding, Zc and 7 can be obtained, which are the fingerprints of the 
winding without knowing the internal detailed physical parameters of the transformer. Zc 
and 7 are further used to detect the winding displacement. 
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The previously shown model is now modified to account for a fault at impedance Z^^ , 
an unknown location in the winding. The winding is split into two sections with effective 
lengths of /, and/^ . The modified model is shown in Figure 8. 

As with the non-faulted winding, the input voltage, input current, and output currents 
5 are known. The two sections on either side of the fault have, in general, different propagation 
constants; for the leftmost section and y2 ^Y^Yx for the rightmost section. Kirchoffs 
voltage law is applied to all of the loops, as was done with the simpler non-faulted winding. 
Four equations result from the four loops in the model. Since there are only four unknowns, 
namely, , // , Zj, , and , the model has a unique solution. The four equations are given 

1 0 by : 

V,-ZJ,=(y^-ZJ,)e-'^ (18) 
(F,+Z<J>"^'=^/+V/ (19) 

iI„R-I„Z,ye-^*'''-^^Vf-ZJj. (20) 

(Vf+ZJ^ =I^^R + ZJ^ (21) 

15 Equations (18) to (21) can be solved numerically to obtam the value of by a 

Matlab script. The corresponding Matlab script for solving equations (18) to (21) is shown in 
Figure 9, It is noted that equations (18) to (21) are defined as functions Fl to F4 respectively 
in the Matlab script. 

Since the per-unit length propagation constant is assumed not to change as the result 
20 of a fault, the location of the fault is simply given by: 

/m = ^-— (24) 

r 

The fault impedance, Z^^ , is also determined by the solution of equations (18) to (21) 
by Matlab. The determination of the fault impedance can give insight into what type of 
25 fault occurred. For example, if is very small, it implies that a short circuit has occurred 
inside the transformer. 
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Anotiier significant parameter is the insulating distance between the windings, called 
internal insulating distance, and the insulating distance between the outer winding and the 
tank, called external insulating distance. The following expressions define capacitance, 
inductance and impedance for a transformer winding: 

C = i^^.lO' F (25) 
4m'Q bd 

Z = 4;r//— JV^^*— .10-* H (26) 

h _b + d ^ ^ 

zc^n^^B^-^ 0 (27) 

WhCTc: 

li\ magnetic permeability of the material 
10 €: dielectric constant of the insulation material 

N: total number of the turns of the winding 
h: axial length of the winding 

b: intemal insulating distance, i.e. distance between the winding and the core 
d:extemal insulating distance, i,e. distance between the winding and the tank 
15 vo: velocity of light in a vacuum, 300m/us. 

As shown by equations (25) and (26), changes in the winding distances (i.e. the 
separation b and d) vary both the capacitance and the inductance. The increase of the 
separation reduces the capacitance between the winding and the ground and raises the self- 
inductance of the winding. Formula (27) further defines that Zc is directly proportional 

20 to . When b is much less than d, which is usually the case for transformers, the 

denominator of the fraction equals 1. Hence the Zc is directly proportional to b. Therefore, the 
intemal insulation distance plays a key role in the value of Zc. When b is in the same order as 
d, botli determine the value of Zc. In conclusion, Zc generally increases with the increase of 
the separation; which can be either the separation between the winding and the core, or the 
25 separation between the winding and the tank. 

It should be noted that expression (27) accurately describes the Zc of large transformer 
windings, since the larger the diameter of the windings, the more accurate formulas (25) to 
(27) become. 
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Therefore, with the above expressions, and after the appropriate measurements are 
conducted^ not only can the characteristic impedance signatm-e value or curve for a particular 
electrical winding be determined, the approximate movement of its windings can be 
determined. 

5 Figure 10 is a graph plotting Zc for an electrical winding with a displacement value of 

b=0, b=0.25 and b=0.5 inches, for a transformer winding alter the application and 
measurement voltage and current across a frequency range of SOOKHz to 2MHz. Curve 100 
represents the Zc curve for a displacement of b=0, curve 102 represents the Zc curve for a 
displacement of b=0.25 inches, and curve 104 represents the Zc curve for a displacement of 

10 b=0.5 inches. As clearly demonstrated, Zc is a smooth function of frequency, and vertical 
shifts of the curve are directly related to the winding displacement. Table 1 below illustrates 
the changes in Zc for radial winding movement in the transformer measurements obtained for 
the plots of Figure 10. It is significant to note that the change in Zc between curves is more 
pronounced at the lower frequency of 500kHz, than at the higher frequency of 2MH2. This is 

15 highly advantageous since the lower frequency signals are less susceptible to noise and other 
electrical phenomenon that can introduce error or inaccuracies into the measurements. 



Table 1 



Frequency 


500 kHz 


2 MHz 


Zc change (for 3% radial 
displacement) 


50 a 


15 S2 


Zc cliaiige (for 6% radial 
displacement) 

* 


75 n 


20 n 



20 Experimental results demonstrate how Zc will change in fault situations, as shown in 

Figure 1 L Figure 1 1 plots five different Zc curves for the same electrical winding, each with 
different fault situations. Curve 110 is the Zc signature for a transformer winding with no 
fault. Curve 112 is the Zc signature when a short circuit happens in the connecting point of 
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the two windings. Curve 114 is flie Zc signature when a fault shunt branch with 220 ohms 
resistance occurs at the same location. Curve 116 is the Zc signature when a fault shunt 
branch with 330 ohms resistance occurs. Curve 118 is the Zc signature when a fault shunt 
branch with 4700 ohms resistance occurs, which is similar to the arcing resistance caused by 
5 partial discharge. It is assumed that the faults are located at the same location as in the short 
circuit case. 

After applying input voltages and measuring output voltages to the winding for each 
case, the Zc results were calculated and plotted in Figure 1 1 . As evidently seen in Figure 1 1, 
in the short circuit case 112, Zc decreased dramatically. Compared with the Zc of the original 

10 curve 110 under normal condition, Zc of the short circuit case 112 is about half of it. In the 
220 and 330 ohms fault cases (curves 114 and 116 respectively), Zc curves cross the original 
curve 110. In the 4700 fault case (curve 118), the Zc curve is above the original curve. It is 
concluded that with the increase of the fault resistance, Zc begins to increase, and crosses the 
original non- fault curve at some certain fault impedances, and continues to go up until it is 

15 above the non-fault curve. Also, short circuit faults can be detected by observing the dramatic 
decrease of Zc. 

Now that the expressions for modeling a transformer winding after a transmission line 
have been described, they can now be applied in methods and apparatus for obtaining a 
characteristic impedance Zc signature from a real transformer. 

20 The characteristic impedance Zc for an electrical winding in a transformer, or 

electrical motor or generator, can be obtained without having to remove the device from its 
native environment. According to an embodiment of the present invention, a self-contained 
portable diagnostic device can be used to obtain Zc of a winding over time. Hence the device 
can be used to detect if there has been winding movement, which could lead to impending 

25 failure of the transformer, or device. 

Figure 12 is a block diagram of a winding diagnostic device according to an 
embodiment of the present invention. Winding diagnostic device 200 is a portable device that 
is used to calculate the Zc of a particular electrical winding based upon the previously 
described equations. Winding diagnostic device 200 includes processing means 202, a signal 

30 generator 204, a sensing means 206, storage means 208, and a user interface 210. Processing 
means 202 is generally responsible for controlling all other components of the winding 
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diagnostic device 200. Those of skill in the art will understaxid that such a device may include 
additional components, which are not shown to simplify the figure. 

Processing means 202 executes algorithms and is progranuned to calculate Zc based 
upon the previously described equations. Processing means 202 can be controlled by user 
5 interface 210 to initiate a measurement function, whereby processing means 202 provides 
output voltage (Vout) and current (lout) information for a specific frequency (f). Preferably, 
processing means 202 steps through a range of frequencies at predetermined step sizes, so that 
a Zc curve can be generated for each specific firequency. Fox each frequency, the sensed 
output voltage (Vout) and current (lout) is used to calculate Zc, and can be subsequently 

10 stored in memory 208. The sensed Vout and lout are preferably received in digital format. 
Once the entire frequency range has been swept, which is preferably between 500kHz and 
2MHz, processing means 202 can plot Zc vs frequency f, to graphically illustrate the signatiure 
of the electrical winding. As will be discussed shortly, processing means 202 can also execute 
Fourier transform algorithms to decompose input and output signals into frequency 

15 components for Zc and gamma calculations for a low voltage impulse test. Ih this test, Vin, 
Tin, Vout, lout are impulse shapes that are recorded and then decomposed into frequency 
components by the Fourier transform algorithm. 

Signal generator 204 can include a function generator, network analyzer, or a simple 
function generator. A network analyzer can produce discrete frequencies one at a time with 

20 high resolution. The simple function generator can inject an analog signal, such as a 

sinusoidal signal at a given frequency. Such signal generators can be used for what will now 
be called a swept frequency test. Alternatively, a recurrent surge generator (RSG) can be used 
to generate low voltage impulse signals for the low voltage impulse test. For example, the 
RSG can generate a train of square shape pulses, or pulses of other suitable shape, which are 

25 injected into the transformer winding. In either case, the ftmction of the signal generator 204 
is to provide an electrical signal having some frequency component. For Zc and gamma 
measurement, the output terminal of signal generator 204 is electrically connected to one 
terminal of the winding. For an on-line transformer, this terminal can be connected to a 
capacitive tap located on the tank of the transformer. Signal generator 204 produces Vout and 

30 lout at a specific frequency fin response to the control signals from processing means 202. 
For off-line applications, current lout that goes into the transformer winding and current lin 
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that comes out of the winding can be measured by measuring the voltage on series resistors. 
For on-line applications currents lout and lin can be measured with current pick up coils 
mounted on the transformer bushings or ottier suitable current detection devices. 

Sensing means 206 can include commercially available high speed digital data 
5 recorders or customized sensing devices. Those of skill in the art will imderstand that such 
customized sensing devices would be constructed with the appropriate resolution for 
converting the analog signals into digital samples that represent the magnitude and phase 
angle of the sensed signals to be used by processing means 202. Those of skill in the art 
understand the functionality provided by commercially available high speed digital data 
10 recorders, such as oscilloscopes for example, and will understand that the core components of 
such boards can be integrated into the presently described diagnostic device 200. Sensing 
means 206 has a terminal connected to the other terminal of the winding, usually another 
capacitive tap of the transformer tank, or a ground bushing terminal, for sensing an analog 
Vin and Tin. 

1 5 Memory 208 can be implemented as any type of known memory, and can be 

embedded with the microprocessor, or implemented as a discrete device. Alternately, memory 
208 can include mass storage media, such as the hard drive of a computer for example. The 
density of memory 208 can be selected to be any suitable size. Memory 208 can be used to 
temporarily store Vout, lout, Vin, In for a number of specific frequencies for subsequent 

20 processing, and/or for storing the calculated Zc for the specific frequency. Those of skill in 
the art will understand that memory 208 can include different types of volatile and non- 
volatile memory, which are addressable using well known methods. 

User interface 210 can include any type of input device and. feedback device, such as a 
screen for displaying alpha-numeric text and graphs. In the present embodiment, user 

25 interface 210 can display a Zc value for a specific frequency, a Zc curve for a range of 
frequencies, and information pertaining to the difference between a baseline Zc value or 
curve, versus a current Zc value or curve. Once the Zc differences have been calibrated to 
corresponding alert thresholds, fiizzy logic can be implemented to display the appropriate 
message for the user. For example, the displayed message can recommend to the user that the 

30 presently measured Zc for a winding indicates potential failure of the transformer, and 
replacement procedures should be initiated. 
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A further discussion of processing means 202 now follows. It is noted at this time that 
processing means 202 can be implemented with commercially available microprocessors, or 
other customizable processing devices such as application specific integrated circuits (ASIC). 
Figure 13 illustrates the preferred functional blocks that are included within processing means 
5 202. 

Processing means 202 includes a basic controller core 250, a sequence generator 252, 
a characteristic impedance calculation engine 254, fiizzy logic engine 256, and a graphic 
engine 258. 

Controller core 250. is responsible for controlling its peripheral components, and their 

10 interaction with each other. Controller core 250 can execute various comparison routines for 
comparing different Zc values or curves to each other and provide a difference value, which 
can then be interpreted by fiizzy logic engine 256. 

Sequence generator 252 can be pre-progranraied with specific Vout, lout values, and 
frequencies f, which are provided to signal generator 204 in accordance with one or more 

1 5 preset sequences. For the low voltage impulse test method, sequence generator 2.52 can 
provide parameters relating to the pulse train to be generated by signal generator 204. 

Calculation engine 254 receives the digital values of Vout^ lout and f, and the digitized 
values of Yin and lin, and is pre-programmed with the transmission line models for a 
winding, as shown by the previously discussed equations and algorithm for executing the 

20 Matlab script shown in Figure 9. The resulting Zc value for a specific input frequency can be 
fed to controller core 250 which can then store the value in memory 208 for later batch 
processing with Zc values for other specific frequencies. Calculation engine 254 can execute 
further mathematical computations, for example, the approximate winding displacement can 
be calculated for a current Zc measurement, relative to a previous Zc measurement. For the 

25 low voltage impulse test method, calculation engine 254 can include the Fourier transform 
algorithms for decomposing applied and received input pulses, where each component 
frequency is then treated as a separate frequency so that Zc and gamma functions can be 
calculated for each one of ttiese frequencies. 

Fuzzy logic block 256 can receive the calculated difference between different Zc 

30 values/curves and apply pre-programmed fuzzy logic to provide corresponding text messages. 
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Alternatively, fuzzy logic block 2S6 can execute the comparison routines instead of controller 
core 250. 

Operating in parallel can be graphic engine 258, which simply plots the calculated Zc 
values against frequency, and presents an exponential Zc signature curve, similar in shape to 
5 those shown in Figure 10. Graphic engine 258 can superimpose more than one Zc curve to 
better illustrate shifts between the curves. 

Since the aforementioned Avinding diagnostic device 200 is portable, any technician 
can quickly and easily obtain a Zo signature of one or more windings of a transformer, or any 
device employing long windings, such as electrical motors or generators. 

10 The cost and complexity of winding diagnostic device 200 can be reduced by 

removing the processing means 202, and consequently, any Zc calculating capability. Hence, 
device 200 would be reduced to a signal generating and electrical parameter 
measuring/collecting device. In such an embodiment, the measured values can be stored in 
memory 208, and downloaded to et laptop computer, or alternatively, sent out through 

15 communication means to a centralized transformer monitoring facility or asset management 
centre. The laptop computer, or the computers at the centralized processing centre, can 
execute all the necessary calculations to generate Zc once the input and output values have 
been provided. In this alternate ennbodiment, the laptop computer, or the computers in the 
centralized processing centre, replace components 202 and 210 of Figure 12, 

20 The previously discussed embodiments employ a signal generator to apply the desired 

input signals to the winding for subsequent output measurements. However, if high freqiaency 
signals are generated within the power system itself during, for example, switching , 
operations, short circuit faults, or lightning discharges, the external signal generator 204 is not 
required. The sensing means can t>e the same as previously described, or alternatively, can be 

25 existing sensors in the system. 

Now a method of calculatixig Zc according to an embodiment of the present invention 
is described with reference to Figure 14. The general method for obtaining the characteristic 
impedance involves applying an input voltage and current to one terminal of the electrical 
winding at a given frequency, and measuring the output voltage and current signal at the other 

30 terminal of the electrical winding. While a single Zc value is sufficient for establishing a 

signature for a particular winding, it is preferable to utilize multiple values to compensate for 
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potential measurement errors. In other words, the use of multiple Zc values, from different 
input voltage/current frequencies, ensures robustness of the overall winding signature. 

The method of Figure 14 assumes that a winding diagnostic device, such as the 
previously described device 200, has its terminals properly connected to the ends of a 
5 winding. The method begins at step 300, where an input signal having a particular Vout and 
lout value is applied to one end of the winding, at an initial frequency f. At step 302, the Vout 
and lout values are measured and stored in memory at step 304. A decision is made at step 
306 to determine if the present frequency is the last frequency of the input signal to be 
applied. If no, then the process proceeds to step 308 where the input frequency is increased by 

10 a predetermined increment of f, and the process repeats at step 300 with the incremented 
frequency. If the present frequency is the last input frequency, then the characteristic 
impedance Zc at each of the input frequencies is calculated at step 310. As previously 
discussed, Zc can be calculated with the previously described equations. Now the 
characteristic impedance signature of tiie winding has been obtained, and can be graphically 

15 plotted to visually verify its exponential shape. The above mentioned method can be modilBed 
to calculate Zc in each frequency loop, and then store each calculated Zc in memory. 

The method of Figure 14 can be modified as follows for implementing the low voltage 
impulse test. At step 300, a train of pulses are applied to the winding. At step 302 the output 
shapes are measured, and then stored at step 304. The method then skips step 306, and 

20 proceeds to step 310 where the stored impiilse shapes are decomposed into frequency 
components, such that the characteristic impedance Zc and gamma functions can be 
calculated for each of the frequencies. Accordingly, steps 306 and 308 are not required since 
the frequency components are obtained through Fourier transform analysis. 

Now that the method for obtaining a signature is described, the signatures can be used 

25 to diagnose the health of the winding. Figure 15 is a flow chart illustrating an example of one 
such method, according to an embodiment of the present invention. 

The method of diagnosing a winding shown in Figure IS starts at step 400, where a 
base characteristic impedance Zc signature is obtained at time ^=0. Preferably, this base Zc 
represents the Zc signature for the winding of a new transformer, which is known to not have 

30 any faults or winding displacement. The base Zc can be obtained through measurement and 
calculation as shown in Figure 14. However, if the transforms has been in operation for a 
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long period of time, a current measurement of its Zc will not be representative of its "new" 
condition. In such a circumstance, it is possible to obtain the base Zc by knowing the physical 
parameters of the transformer winding, and applying equation (27) above. Once obtained, the 
base Zc is stored in memory for that particular transformer winding at step 402. At a 
5 subsequent period in time, another measurement of Zc is made at step 404. This measurement 
can be done at regular intervals, such as once every year, or shortly after a thunderstorm. This 
current Zc is stored in memory for the same transformer winding at step 406. Having the 
current Zc and the base Zc, a difference between the two can be easily calculated at step 408. 
The magnitude of the difference can indicate the level of winding movement, where a 

10 relatively large difference can indicate potential for transformer failure, while a negligible 
difference indicates very little to no winding movement. Of course, a complete history of Zc 
for a particular winding can be plotted to show progressive winding movement over time. An 
extrapolation of the Zc history may indicate the failure potential of the transformer in the 
future. Accordingly, the transformer can be replaced before catastrophic and dangerous failure 

15 occurs. 

An alternate technique for diagnosing the relative health of a winding can be 
conducted as follows. Standard three phase transformers have three identical windings. Since 
it is unlikely that all three windings will suffer identical movement, the Zc for each winding 
can be measured and calculated, and then compared to each other. If the Zc of all three 

20 windings are substantially similar, then there has not been any significant movement in any of 
the windings. However, if one Zc is different than the other two, then the technician can 
quickly identify the winding that has suffered some type of winding movement. 

While the previously described embodiments of the present invention illustrate a low 
input signal frequency of 500kHz, even lower input signal frequencies can be used to 

25 calculate the characteristic impedance of the electrical winding. The minimum frequency 
input signal that can be used for any given electrical winding will be dependent on the 
maximum length of the winding, which can easily be determined by any person of skill in the 
art. As previously mentioned, if a sufficient number of wavelengths of the signal can 
propagate through winding, then the wave propagation functions can be used to describe it. 

30 Naturally, the data recorders should preferably be highly accurate and relatively immune to 
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electrical noise and other phenomenon that can adversely affect the output signal 
measurements. 

While the previously described embodiments of the present invention describe 
methods and apparatus for obtaining characteristic impedances of electrical windings of 
5 transformers that are cn-hne^ meaning that they are connected to the main electrical system, 
they can also be applied to off-line transformers. A disconnected transformer being tested 
according to the embodiments of the invention will yield characteristic impedance values 
which are at least as accurate as those determined while the transformer is on line. 

By taking advantage of the wave propagation property and frequency-dependent 
10 transmission line model, one can simulate the real transformer winding as a high frequency 
distributed transmission line, obtain the characteristic impedance and propagation constant of 
the winding, calculate the travelling time for certain frequency for the total equivalent length 
of the winding. By comparing the historical change of the characteristic impedance Zc the 
presence of winding distortion or movement can be detected. The fault location results in an 
15 equivalent sub length, which can be compared with the total equivalent length and ftirther 
used to determine where the winding movement or distortion has occurred. 

The above-described embodiments of the present invention are intended to be 
examples only. Alterations, modifications and variations may be effected to the particular 
embodiments by those of skill in the art without departing from the scope of the invention, 
20 which is defined solely by the claims appended hereto. 
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